Technique for detecting warm-hot intergalactic gas 
in quasar UV spectra 
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ABSTRACT 

The ionizing spectral energy distribution of quasars exhibits a steepening of the distribution 
shortward of ^ 1200 A. The change of the power-law index from approximately —1 (near-UV) to 
—2 (far-UV) has so far been interpreted as being intrinsic to quasars. We consider the possibility 
that the steepening may result from a tenuous absorption component that is anticorrelated with 
large mass overdensities. UV sensitive satellites, whose detectors can extend down to 1000 A, can 
set a useful limit to such an absorption component through the search of a flux increase in the 
window 1050-1190 A (observer frame) with respect to an extrapolation of the continuum above 
1230 A. Since the recent FUSE or HST-STIS data do not show any obvious discontinuity in this 
region, this effectively rules out the possibility that intergalactic Hi absorption is very important, 
and it is concluded that most if not all of the steepening is intrinsic to quasars. A smaller flux 
discontinuity of order 1% cannot, however, be ruled out yet and would still be consistent with 
the warm-hot intergalactic component if it amounts to 30% of the baryonic mass, as predicted 
by some models of large scale structure formation, provided its temperature lies around 10^ '^ K. 



Subject headings: galaxies: intergalactic medium ■ 
— radiative transfer — ultraviolet: general 

1. Introduction 

The ionizing spectral energy distribution (here- 
after ISED) of nearby active galactic nuclei can- 
not be observed directly, due to the large Galactic 
absorption beyond the Lyman limit. Owing to 
the redshift effect, however, we can get a glimpse 
of the ISED from the spectra of distant quasars. 
The pioneering work of Zheng et al. (1997, ZK97), 
using HST-FOS archived data, identified a notice- 
able steepening in quasar SEDs near 1200 A. The 
power-law index a (F^, (x i>°') in quasars steepens 
from —1 for A > 1200 A to ~ —2 at shorter wave- 
lengths^, down to 350 A (Xrest)- In a more recent 
study. Teller et al. 2002 (hereafter TZ02) found 
similar results with a composite spectrum charac- 
terized by a mean near-UV index of —0.7 steepen- 



^Throughout the text, A^jj and \reat will denote wave- 
lengths in the quasar rest-frame or in the observer-frame, 
respectively [Xrest = (1 + zq)~''- \of,s.]- 



large-scale structure of universe — galaxies: active 

ing to ~ —1.7 in the far-UV. The above authors 
favor the interpretation that the steepening is in- 
trinsic to quasars and that it is the signature of a 
comptonized accretion disk. 

On the other hand, certain distributions of in- 
tergalactic absorption gas can also cause an (ap- 
parent) steepening of the SED, starting at 1200 A 
(Arest), as shown in this Paper, that is, located 
in the same position as the break encountered by 
TZ02. As long as ISED studies are based on de- 
tectors that do not extend beyond 1200 A (Xobs.), 
one cannot readily distinguish between the con- 
tribution to the observed steepening from an in- 
tervening absorption model and that of a purely 
intrinsic break in the quasar SED as proposed by 
ZK97 or TZ02. FUSE or HST-STIS spectra, how- 
ever, extend much farther into the UV and will be 
shown here to provide us with a compelling test 
to discriminate between the two interpretations. 

Large scale structure formation models predict 
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that a substantial fraction of baryons should reside 
in a warm-hot phase at current epoch, possibly up 
to 30-40% of Ofear. (Dave et al. 2001). Depending 
on its temperature and distribution with redshift, 
this warm-hot gas may contribute to a reduced 
fraction of the observed steepening of the ISED. 
In this Paper, we show that this component may 
give rise to a slight flux increase (that is, a dis- 
continuity) in the region 1050-1 190 A (Xobs.) with 
respect to longer wavelengths. The calibration of 
the level of this discontinuous excess flux can be 
used as a technique to set useful limits for the 
baryonic mass contribution from the warm-hot in- 
tergalactic component. 

The objective (and results) of this Paper is 
threefold: 

1. Find the simplest and yet physically mean- 
ingful absorption gas distribution which can 
mimic the observed steepening (§3.1), 

2. Devise a technique to falsify the previous 
proposition that the break is due to absorp- 
tion, thereby confirming the idea that the 
steepening is intrinsic to quasars (§3.2), 

3. Propose a dependable technique to reveal 
traces of a warm- hot baryonic component in 
the local Universe using UV spectra (§3.3). 

2. Procedure and calculations 

2.1. Generalized Gunn-Peterson effect 

The Gunn-Peterson (GP) effect (Gunn & Pe- 
terson 1965) sets stringent limits for the presence 
of neutral diffuse gas at high redshifts. In the 
simplest form of the GP test, the absorption gas 
produces a flux decrement between Lya and Ly/3. 
The decrement is measured against a continuum 
level, usually taken to be a power-law extrapo- 
lated from the region redward of the Lya emission 
line. With high resolution spectra, this technique 
is sensitive to H l columns otherwise too small to 
produce resolvable Lya lines. Songaila, Hu, & 
Cowie (1999) in this manner could set a limit of 
the GP opacity tgp < 0.1 at z = 4.7. Another 
manifestation of the GP effect is illustrated by the 
staircase transmission curve of M0llcr & Jakobsen 
(1990), which is the result of unresolved Lya for- 
est lines in intermediate resolution spectra. This 
transmission curve is also characterized by a broad 



trough (partly due to photoelectric absorption), 
the Lyman valley, whose minimum occurs at a 
rest-frame wavelength ~ 650 A. Before averaging 
their quasar spectra, ZK97 and TZ02 statistically 
corrected each quasar spectrum for the presence 
of such a Lyman valley by calculating the appro- 
priate transmission curve, using the scheme devel- 
oped by M0ller & Jakobsen (1990). Finally, the 
GP effect can be generalized to include the broad- 
ened wings of faint Lya forest lines as a result of 
the Hubble flow in underdense regions. 

Our contention is that these manifestations of 
the GP effect are sensitive to certain classes of H i 
distributions as a function of redshift. Other Hi 
distributions might instead be associated to large 
scale voids. In that case, the amount of absorp- 
tion gas would be decreasing towards the back- 
ground quasar, against which we are trying to de- 
tect it. Such opacity behavior might be harder to 
pinpoint, since it would not necessarily produce 
a recognizable absorption break at the rest-frame 
Lya. Another factor that could hinder the detec- 
tion of intergalactic absorption is that the absorp- 
tion component evolves in the opposite direction 
as the Lya clouds, becoming in other words pro- 
gressively more abundant towards lower redshifts. 
For instance, hierarchical structure formation is 
expected to generate a hot gas phase, whose mass 
is predicted to increase with time (e.g. Phillips, 
Ostriker, & Cen 2001), that is, with a temporal 
evolution which is opposite to that of the Lya for- 
est. Any residual absorption it may cause would 
occur at much shorter wavelengths, in a wave- 
length domain too far from the Lya emission line 
(Arest) to rely on an extrapolation of the contin- 
uum observed redward of this line. 

With these considerations in mind, we have ex- 
plored various H l distributions, whose absorption 
behavior is not as evident as that expected from 
the traditional GP effect. The hypothesis being 
tested is this Paper is whether all or just a small 
fraction of the far-UV steepening in Fi, is caused 
by absorption. We will demonstrate in § 3 how 
this hypothesis is observationally falsifiable. Even 
if most of the SED break turns out to be intrinsic 
to quasars, the technique can still be used to set 
a lower limit on the temperature of the warm-hot 
intergalactic medium (WHIM). 
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2.2. Derivation of the transmission curve 

The technique used to genc;rate a simulated 
composite spectrum to compare with the observed 
spectrum of ZK97 is described in Binette et al. 
(2002). Briefly, wc assume a typical spectrograph 
sensitivity window for FOS of 1300-3000 A. We 
then simply redshift this window in locked steps 
in order to simulate quasars of different redshifts 
in the range 0.33 < zq < 3.6 (the same range as 
ZK97). Before averaging, each redshifted quasar 
SED is multiplied by the redshift integrated trans- 
mission curve Tx (defined below). We assume the 
concordance ACDM cosmology with Sl\ = 0.7, 
f^M = 0.3 and h = 0.67 with h = Fq/IOO. 

To be definite, we adopt a uniform gas distri- 
bution (see Binette et al. 2003). Equivalent cal- 
culations assuming a clumpy medium were also 
carried out for comparison. We found that the 
ill-defined additional parameters required then do 
not reveal anything fundamental or even interest- 
ing. Furthermore, the two formalisms are equiva- 
lent for the assumed opacity regime in which the 
diffuse absorption component lies in the linear part 
of the curve of growth, that is for Hi columns 
Nho < 10^2 (.ni-2. 

The simulated quasar spectra are divided in en- 
ergy bins, and for each quasar rest-frame wave- 
length bin j, we calculate the transmitted inten- 
sity P{, = hjTxj = h^e-'^^^^^ by integrating the 
opacity along the line-of-sight to quasar redshift 

T(A,) = El=oC^i(iqi)«Ho(^)f 

where is the quasar rest-frame wavelength for 
bin j and Uho (z) the intergalactic neutral hydro- 
gen density, which consists of one of the three 
density distributions discussed below [eqs (1—3)]. 
The summation is carried out over the follow- 
ing opacity sources: photoionization (i = 0) and 
line absorption from the Lyman series of hydrogen 
(1 < i < 10). Although our code could include up 
to 40 levels, we found that considering only the 10 
lowest proved to be adequate. We adopted a fidu- 
cial velocity dispersion b of 30 km/s and assumed 
a Gaussian profile for the cTj of the lines. 



2.3. Heuristic H i geis distributions 

2.3.1. Three Hi distributions to test 

In calculating the transmission function of a 
quasar at redshift zq, we will consider three dif- 
ferent distributions of absorption gas density, uho , 
with redshift, which we will then compare with the 
composite SED of either ZK97 or TZ02. The first 
is given by 

nHo{z) = n%o{l + zf{l + zr{z'Q)!' (1) 

where z is the absorbing gas redshift, zq the 
quasar redshift, z'q the quasar redshift as seen 
from the absorbing gas at z [that is = (1 -|- 
zq)/{1 + z) — 1], and n^o the neutral gas density 
at zero redshift. 7 and /3 are adjustable parame- 
ters. The {1 + z)^ factor represents the cosmolog- 
ical expansion of the Universe. The factor {z'q)^ 
with /? > implies that the density decreases in 
the neighborhood of individual quasars. 

The second distribution which will be consid- 
ered is 

nHo{z)=n%o (l + z)^ (2) 

where zp is the size of the region (near each 
quasar) within which the density decreases as 
{zp/z'qY towards the quasar [note that /3 > 
as in eq. (1)]. It is possible to substitute the de- 
nominator above by the expression 1 + [rp/r'q)^ . 
For each value of zp discussed below, we we have 
searched for the value of rp which would give us an 
equivalent fit of the composite ISED. The param- 
eter Tp represents the zone of influence (or cav- 
ity radius) of the quasars in spatial units. Al- 
though the density decreases smoothly towards 
each quasar, in this text we will refer to this zone, 
where r'q < rp near each quasar, as a 'cavity'. 
The third distribution is given by 

n^o{z)=n%.{l + zf^-%^^}^ (3) 

where the exponential function in the numerator is 
a parametric fit to the increase of the warm-hot in- 
tergalactic medium (WHIM) towards the current 
epoch as calculated by Dave et al. (2001) (their 
Model D2). We will again discuss which value of 
rp produces a fit equivalent to that provided by 
our favored zp values. 
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2.3.2. Three conditions to satisfy 

The above distributions lie in a sequence of in- 
creasing plausibility, as discussed below. We em- 
phasize that our exploration included many more 
functions than the above three, the vast majority 
of which could not reproduce a smooth steepening 
of the spectrum near 1100 A (Xobs.) (failing condi- 
tion a below). Nevertheless, these three distribu- 
tions will suffice for the purpose of unequivocally 
testing the absorption hypothesis. 

In our search, the following procedure was 
adopted: we explored a large set of functions with 
the following three priorities in mind: a) the Uho 
distribution must produce an acceptable fit of 
the the composite ISED without a discontinuity 
near Lya, b) it must be a physically meaningful 
function, and c) it must be consistent with other 
information we have about the Universe such as, 
for instance, the volume density of quasars. We 
were able to satisfy conditions a and b, although 
c could not be strictly respected for our favored 
test case Model C, as discussed below. 

3. Model results 

3.1. Reproducing the observed steepening 

If we adopt the first distribution of H i corre- 
sponding to eq. (1) and the parameters listed in 
Table 1 for Model A. that is, a zero redshift density 
n%o = 4.7 X 10-12 cm-'\ 7 = -1.5 and /? = 0.8, 
we can reproduce the observed smooth steepen- 
ing of the composite spectrum of ZK97 remark- 
ably well as shown in Fig. 1 by the thick solid 
line. The intrinsic SED considered in Model A is 
a single power-law of index —1 (short-long dashed 
line), which we extend to all wavelengths. This in- 
dex provides a good fit at the longer wavelengths 
(1300-2300A) of the underlying continuum in the 
ZK97 composite spectrum. 

This first distribution has also the property that 
the far-UV index, apuv, defined in the range 350- 
900 A (Xrest ), increases very little with quasar red- 
shift (Binette et al. 2003). This agrees with the 
findings of TF02 (c.f. their Fig. 12). However, al- 
though the simplest, this distribution is not phys- 
ically meaningful in the following aspect: the lim- 
itless range of the dependence of eq. (1) on z'q 
implies that the local Hi density depends on the 
quasar redshift. To clarify this point: two inde- 



pendent quasars, one nearby and the other at red- 
shift 3, but lying along very nearby line-of-sights 
on the sky, would be characterized by very differ- 
ent local Hi densities according to eq. (1). Such 
behavior is unphysical, because the Hi density in 
the local Universe should not depend at what dis- 
tance the background quasar is located. We there- 
fore conclude that condition b ('physically mean- 
ingful', see §2.3.2) is not satisfied and that the 
rijjo distribution of eq. (1) must be rejected. In 
the text, we will refer to this undesirable property 
as the "nearby line-of-sight problem" . 

A distribution without any dependence on z'q 
would not suffer from this problem. However, 
in all the distributions that succeeded in repro- 
ducing the smooth roll-over near 1200 A, some 
kind of smooth gas cavity near each quasar had to 
be taken in, otherwise, a deep and sharp absorp- 
tion discontinuity^ unavoidably appears at 1216 A 
(Arest) in the synthesized composite SED, a fea- 
ture that is obviously not observed. The thin long- 
dashed line in Fig. 1, for instance, illustrates the 
magnitude of this sharp discontinuity present near 
Lya when no cavity is considered (here with (3 set 
to zero in eq. (1)). 

The solution to this problem is to have the de- 
pendence on z'q limited to a reduced zone of influ- 
ence near each quasar. This is achieved with the 
expression in the denominator in eqs (2) and (3), 
where zp limits the distance scale within which 
the density is decreasing (towards the quasar) as 
a power-law of index f3. We therefore have two 
competing effects present in our subsequent dis- 
tributions [of eq. (2) or eq. (3)]: one factor de- 
scribes the H I density as a function of redshift 
(seen from us) z, which has no other QSO close to 
the selected line-of-sight, and the other describes 
the H I density within the neighborhood of QSOs 
(within a zone of influence c^i rp). Two differ- 
ent line-of-sights lying at a projected distance on 

^In Binette et al. (2002), such a discontinuity at 1216 A 
(Areat) was Smoothed out by invoking photoionization by 
the background quasar (in analogy with the "proximity ef- 
fect"). However, the metagalactic background radiation 
(relative to a single quasar) is much too strong for this 
proposition to be sustained. Eastman, MacAlpine, & Rieh- 
stone (1983) avoided the discontinuity by invoking absorp- 
tion cloudlets being progressiveh' accelerated up to 0.8c 
within a region < 10 kpc near the background quasar. The 
main problem was that the cloudlets themselves emitted 
more UV than the background quasar. 
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the sky greater than rp will be described by a 
unique and self-consistent distribution of H i with 
redshift. As a consistency check, the volume occu- 
pied by all the QSO cavities must be smaller than 
the total volume sampled. 

A possible interpretation of the decreasing den- 
sity towards individual quasars is that our puta- 
tive Hi distribution is associated to the largest 
scale voids. The environment of quasars might 
also be hotter and transparent as a result of (pro- 
tocluster) stellar winds from the associated large 
scale and most massive structures. Our aim will be 
to make zp as small as possible, but without pro- 
ducing any discontinuity or dip near 1200 A. The 
short-dashed line in Fig. 1 shows that the distri- 
bution corresponding to eq. (2) (Model B) results 
in as good a fit to the ZK97 spectrum as Model A, 
when using 7 = —3.0, zp = 0.7 and (3 = 1.5 (see 
Table 1). 

Nevertheless, this Model B is not satisfactory. 
In effect, a cavity with zp = 0.7 requires a cavity 
of size rp ~ 2 Gpc, which is more than half the 
path-length to the highest-z quasar in the sam- 
ple (i.e. 3.75 Gpc for zq = 3.6). Hence most 
line-of-sights would cross many overlapping cavi- 
ties (hence condition h is not satisfied) and such a 
situation not properly described by eq. (2). Con- 
dition c ('consistency with other reliable informa- 
tion', see §2.3.2) is not respected either. First, 
the cumulative volume of quasar influence on H i 
exceecds 15% of the total volume sampled. Sec- 
ond, the exponent of (1 + z). 7 = —3.0, is much 
too steep if we were to compare our Model B to 
model-predictions of the evolution of the WHIM. 
For instance, most models in Dave et al. (2001) 
are characterized by redshift-averaged values of 7 
in the range —1.1 to —1.5. 

In order to ensure a physically more meaningful 
H I distribution, we replaced the (7) power-law de- 
pendence on redshift by an exponential fit to the 
WHIM Model D2 of Dave et al. (2001), the de- 
scription of which lies in the numerator of eq. (3). 
Replacing the previous best fit distribution with 
7 = —3 by a distribution that matches the behav- 
ior of the WHIM with redshift better [by using 
either eq. (3) or 7 ~ —1.2 and eq. (2)] comes at a 
price, however. In effect, such distributions always 
result in a smooth rise in transmission towards the 
short wavelength extremity of the SED, beyond 
500 A. Interestingly, the more recent composite 



spectrum from TZ02 presents such a rise in the far 
UV as displayed in Fig. 2. TZ02 have questioned 
the reality of this rise. The latter might be caused 
by a tendency of high redshift quasars to have a 
harder ISED. Nonetheless, because the TZ02 work 
contains more objects than the ZK97, and because 
this feature is seen in both radio-loud as well as 
radio-quiet quasars, we will take their results at 
face values and compare their radio-quiet com- 
posite SED with models which used eq. (3). The 
thick solid line in Fig. 2 represents our Model C, 
which is characterized by a significantly smaller 
zp = 0.3. Although the fit is imperfect, it matches 
the overall trends present in the TZ02 composite 
spectrum. The adopted intrinsic SED used here 
is somewhat harder, corresponding to an index of 
-0.72 (which is the value inferred by TZ02). 

For Model C, with zp = 0.3 we find that 
rp ~ 800 Mpc. A single cavity occupies there- 
fore 1% of the total volume sampled. Using the 
known luminosity function of bright QSOs [i.e. a 
quasar density of order ^ = SxlO"*" Mpc~^ (Boyle 
2001)], the volume"^ corresponding to this value of 
rp is too large by a factor ~ 650 and condition c is 
therefore not fulfilled on the ground of the known 
density of QSOs. With this caveat in mind, we 
proceed to study the behavior of the absorption 
of our putative intergalactic component at wave- 
lengths shorter than provided by HST-FOS. This 
will allow us to calibrate a new method to detect 
the WHIM and, indirectly, to determine to what 
extent condition b is satisfied in Model C. 

3.2. Prediction of a possible discontinuity 

in the far UV 

3.2.1. Analysis of the transmission curves 

The derivation of the composite SED of Model 
C required using eq. (3) and calculating the trans- 
mission function at evenly spaced redshift values. 
We now turn to the particularities of these trans- 
mission curves. Two such transmission curves are 
presented in Fig. 3 for illustrative purposes for 
quasars of redshifts 2.76 and 0.33. Clearly, the 
main opacity source is the Lya line. The con- 
tribution of photoelectric absorption to the total 
opacity is relatively small as illustrated in the Fig- 

*That is, one should realize the condition Vq = |7rrp^ < 
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lire while the contribution of higher series hnes can 
be appreciated by the reduced size of the jumps 
that are visible in the transmission curves. The 
main result to emphasize is that a discontinuity 
appears to the blue of Lya {Xobs. ) (that is beyond 
the UV limit of the HST-FOS detector) in all the 
calculated transmission curves, not only when as- 
suming the parameters of Model C, but in all the 
Hi distributions that satisfied condition a ('good 
fit', see §2.3.2). 

By comparing the transmission curves in Fig. 3 
with those resulting from Lya forest absorbers, a 
clear difference emerges. The transmission func- 
tions derived by (M0ller & Jakobsen 1990) (or 
ZK97) are characterized by a stair-case behavior, 
in which the opacity from each physical process 
(either a Lyman series line or photoelectric ab- 
sorption) decreases towards shorter wavelengths, 
due to the strong evolutionary nature of the Lya 
forest, in which the density of absorbers increases 
steeply with redshift. In the case of the transmis- 
sion resulting from equations (1) (3), the opacity 
(along a given step of the staircase transmission 
curve) tends to either remain high or increase to- 
wards shorter wavelengths (but, of course, only up 
to the threshold wavelength of the absorption pro- 
cess involved). The obvious reason for this marked 
differenc-c! is that the H i distributions used here are 
characterized by a density that increases towards 
smaller redshifts z. 

An inspection of Fig. 3 suggests that, rather 
than looking for a sharp absorption edge near Lya 
in the quasar rest-frame, the most obvious demar- 
cation produced by n^o distributions (which sat- 
isfy condition a of a 'good fit') is at the shorter 
wavelength end, in the 1190 1070 A observer- 
frame region where Ly/3 rather than Lya domi- 
nates the absorption. In order to detect such a 
discontinuity near 1216 A (Xobs.), one requires de- 
tectors with a sensitivity, which extends to wave- 
lengths shorter than the earlier HST-FOS window, 
such as provided by FUSE or HST-STIS. 

3.2.2. Operative measurement of the jump: tuqq 

In order to test our absorption hypothesis 
against individual quasar spectra observed by 
STIS or FUSE, we must first calibrate the depth 
of the expected discontinuity (on the blue side of 
local Lya) as a function of n^o for our different 
models. Since the intrinsic far-UV spectral index. 



c^FUVi is not precisely known and, as suggested 
by the substantial scatter found by TZ02, may 
even vary from quasar to quasar, it makes sense 
to measure the depth of the jump in a way that 
does not depend on an a priori knowledge of the 
intrinsic ISED. The technique proposed here is 
to evaluate the discontinuity's depth at 1160 A 
{Xobs.) by comparing the flux there to the extrap- 
olated value from a power-law fit redward of Lya 
(Xobs.), within the narrow window 1260 1360 A. 
This window is meant to exclude the Galactic Lya 
absorption trough as well as the geocoronal Lya. 
We define the quantity tubo = logo [F^bs" / ^exfr.]^ 
which can be shown to be insensitive to the power- 
law index of the intrinsic ISED assumed. In Fig. 3, 
we illustrate the procedure by taking as example 
the transmitted SED of a. zq = 1.0 quasar and 
Model C (the calculation assumes an ideal detec- 
tor without wavelength coverage limitations). The 
two squares illustrate the position (albeit here in 
the quasar rest-frame), at which it is proposed 
to define Tugo. For each model discussed (A-D), 
we list in Table 1 the value of Tugo evaluated at 
Zq = 1, while in Fig. 4 we plot the behavior of 
T1160 as a function of zq. 

3.2.3. Condition b and the behavior o/rneo 

As discussed in §3.1 concerning the "nearby 

line-of-sight problem" , it is unphysical to have the 
local density vary when the background quasar, 
for instance, lies at redshift 2 rather than 3. We 
can easily find out when a distribution suffers from 
this problem by checking whether or not rneo is 
constant at moderate and high values of quasar 
redshift. It is apparent in Fig. 4 that Models A and 
B suffer severely from the above-mentioned prob- 
lem. In the case of Model C, however. Turn is flat, 
at least beyond zq > 0.6. However, within a ra- 
dius around us given by the large value of rp ~ 800 
Mpc, we may reasonably expect to lie within the 
zone of influence of a (dominant) quasar. In this 
case, rather than the smoothly increasing function 
depicted in Fig. 4, a single value of rffgo may in- 
stead apply (riigo reflects the local Hi density). 
Its value would depend on our distance from this 
dominant quasar. Furthermore, nieo would not 
necessarily be isotropic. In conclusion, the smooth 
initial rise of rneo in Model C is at best an ideal- 
ization. Given the likelihood of being positioned 
inside the cavity of a single quasar, the predicted 
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Tiieo constitutes an upper limit (a more probable 
value may lie a factor of a few lower). 

3.2.4- Negative results using two EUV quasar 

spectra 

We now use Model C as a test case, since it 
is the most physical at hand under the assump- 
tion that the 500 A rise in the composite ISED of 
TZ02 is real and that the break is entirely due to 
intergalactic absorption. For any quasar zq > 0.5, 
one finds that Tugo is quite significant. Typically, 
Tiieo > 0.25 as seen in Fig. 4. If that much absorp- 
tion were present, it would therefore be a striking 
feature in the far-UV spectra [provided the quasar 
could be observed blueward of 1200 A {Xobs )]■ It 
turns out, however, that the observations of the 
quasar HS1543+5921 by (Bowcn, Tripp, & Jenkins 
2001, ZQ = 0.807) with STIS and HE2347-4342 by 
(Kriss et al. 2001, zq = 2.885) with FUSE, show 
no indication at all of a discontinuity at the ex- 
pected level of T1160 = 0.27 and 0.28, respectively, 
m the 1160 A region, as compared to the extrap- 
olation from the 1270 A region. We consider that 
a discontinuity as small as rneo = 0.05, if present 
in the two spectra, would have been apparent. 

We (;onchidc that intergalactic absorption can- 
not be the main cause of the observed steepen- 
ing in the composite SED. Model C fails not only 
because of the observed absence of a discontinu- 
ity in the far-UV but also because it is inconsis- 
tent with the known density of quasars (condition 
c) . No other H l distribution could be found that 
would solve alll three conditions. Our study there- 
fore confirms that the change of slope near 1100 
A must be in origin intrinsic to quasars, as pro- 
posed by ZK97 and TZ02. This general conclu- 
sion is not dependent on the particular distribu- 
tion adopted, since in all the distributions we have 
explored (which satisfied condition a), Tneo was 
always > 0.2 when zq > 0.5. 

3.3. Detecting the WHIM in absorption 

Having recognized that most of the break is in- 
trinsic to quasars, we now turn to the problem 
of determining how much flux discontinuity near 
1160 A (i.e. Tiieo) can be expected if we assume 
the H II density predicted by WHIM models and a 
fiducial temperature of lO^ "^ K. For this purpose, 
we adopt an intrinsic SED, which already incorpo- 



rates the intrinsic far-UV steepening. It consists 
of a broken power-law, which has the same index 
of —0.72 in the near-UV as in Fig. 2, but that 
sharply turns over at 1200 A (Xrest) into a steeper 
index of —1.57 in the far-UV, as represented by 
the short-long dashed line in Fig. 5. These two in- 
dices correspond to the values characterizing the 
radio-quiet quasar composite of TZ02 while the 
turn-over wavelength is within the range of values 
encountered by TZ02 (w 1200-1300 A). 

3.3.1. Model D using a broken power-law SED 

Adopting eq. (3) and the above broken power- 
law, we first determine which values of n^o and zp 
best reproduce the observed composite SED. The 
result is Model D, plotted as the solid line in Fig. 5. 
The same /? = 1.5 is used as for Model C and zp 
could now be set to the small value of 0.08 with- 
out producing any significant discontinuity* on the 
blue side of the quasar Lya emission line. Model 
D succeeds rather well in fitting the TZ02 com- 
posite. It does not constitute, however, a unique 
solution because of the uncertainties regarding the 
intrinsic SED. If for instance, we shifted the break 
from 1200 to 1100 A, the assumed ISED would lie 
higher above the composite of TZ02 and a higher 
density WHIM would therefore be required for the 
model to overlap the data. 

With a predicted value of Tneo as small as 0.05 
(see Table 1), Model D is characterized by a much 
smaller discontinuity than Models A-C, simply be- 
cause the assumed ISED is much closer to the 
observed composite of TZ02. Furthermore, rneo 
does not depend on redshift anymore, as shown in 
Fig. 4, which indicates that the cavity near the 
individual quasars is not affecting our determina- 
tion of n^o ■ Inspection of the published spectra of 
HS15434-5921 by (Bowcn, Tripp, & Jenkins 2001) 
and HE2347-4342 by (Kriss et al. 2001) suggests, 
however, that rfflg is likely to be smaller than 
0.05. A thorough analysis aiming at setting strin- 
gent upper limits on rf^lo would clarify the matter. 

*The value of rp which can provide an equivalent fit is 200 
Mpc. Interestingly, this value is only a factor of two too 
large with respect to the density of quasars at high red- 
shifts. This can be resolved if we use a function rp of red- 
shift such as 300 (l-t-z)"! Mpc. Although r p is then larger 
in the local Universe, at redshifts ~ 2 we get rp ~ 100 
Mpc, which is the epoch when the density of quasar is 
known to peak. Condition c is marginally satisfied in Model 
D. 
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3.3.2. Deriving rneo from the expected WHIM 

Let us now turn to the inverse problem of deter- 
mining the height of the jump rneo expected from 
a WHIM component that would contribute 30% of 
the total baryonic mass in the local Universe (as 
in Model D2 of Dave et al. 2001). One advan- 
tage of the operative definition of rneo proposed 
in §3.2.2 is that the choice of a given intrinsic SED 
is not critical to the purpose of determining n^o- 
This is verified from Table 1 where the quotient 
n^jo/Tiwo evaluated at zg = 1 in Models C and D 
are 1.89 x 10"" and 1.80 x 10"" cm'^, respec- 
tively. They diff'er by only 5%, despite the use of 
very different intrinsic SEDs. Because n^o scales 
linearly with Tngo, all the information needed to 
derive rneo can be extracted from Model D inde- 
pendently of whether or not the density inferred 
turns out very different than in Table 1 or whether 
the ISED assumed is somewhat incorrect. 

For any arbitrarily small observed discontinuity 
rfilo, the inferred Hi density is given by n^jo = 
1.8 X 10~^^ '''iilo cm~^ which does not depend on 
zq because we have zp ^ 0.4. In terms of the crit- 
ical density SHQ/SnGniH, this density becomes 

O _ 1.6x10"" _obs. (A\ 

^^WHIM — — j-p — T^iieo 

where x is the neutral fraction of the putative 
WHIM absorption gas component. We can ex- 
press ^wHiM as a fraction Fb of the (current 
epoch) total baryonic mass 0,bar. = 0.021/i~^ 

Fb = nwHIM/flbar. = '^^^ T^rib (5) 

With Fb = 0.3, we obtain that T^^flo ~ x/2.5 x 
10"^. Since the WHIM is Hkely to be highly 
ionized, we are lead to expect very small val- 
ues of T°iQQ. For example, for a collisionally ion- 
ized WHIM at a temperature 10^ '^ K, we derive 
a; ~ 3.1 X 10-6 and therefore rfflo = 0.012. 
Clearly, a high S/N spectrum is required in order 
to be able to detect a 1% excess flux at 1160 A. 
It remains nevertheless a worthwhile exercise to 
carry out, since any limit that can be put on t°iqq 
would directly translate into a lower limit on the 
neutral fraction (hence on the temperature) of the 
WHIM component. Note that, with that small 
TiiQQ, we have the liberty to set zp to zero in eq. 
(3) and condition c (§2.3.2) is therefore implicitly 
satisfied. By the same token, we do not have to 
worry about whether we possibly lie in the zone 



of influence of a single quasar anymore. In order 
to optimize the detection sensitivity, the strategy 
to follow would be to add together as many EUV 
quasar spectra as possible (ideally with redshift in 
the range 0.5 ^ zq < 0.8) and then extract rfflo 
from the summed stack. 

4. Conclusions 

In all the distributions of Hho {z) used for re- 
producing the steepening observed in the energy 
distribution of quasars, the dominant atomic pro- 
cess contributing to the opacity is Lya scattering. 
An important implication of these calculations is 
that the transmission curve shows a significant dis- 
continuity to the blue of Lya {Xobs.) where the 
transmission rises sharply. The size of the discon- 
tinuity, T1160, for any distribution of intergalactic 
gas that can successfully reproduce the SED steep- 
ening reported by TZ02 (assuming a single power- 
law SED), is predicted to be typically > 0.25 in 
spectra of quasar redshift zq > 0.5. Since the two 
quasars HS1543+5921 and HE2347-4342, which 
were observed by HST-STIS and FUSE, respec- 
tively, do not show any sign of a discontinuity any- 
where near this level, we consider the case proven 
that most if not all of the 1200 A continuum steep- 
ening is intrinsic to quasars rather than due to H i 
absorption. On the other hand, a small fiux in- 
crease by as little as 1% percent is still possible 
and would be consistent with a WHIM tempera- 
ture ~ 10^'^ K contributing up to 30% of the cur- 
rent baryonic matter. 
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Table 1 
Model parameters 



Label 


Equation 
(#) 


sed 

A < 1200 A 


<o(2: = 0) 
(cm-3) 


7 





zp 


T1160 
(^Q = 1) 


A 

B 
C 
D 


(1) 
(2) 
(3) 
(3) 


^-1.57 


4.7 X 10-1^ 

1.5 X 10"" 
5.2 X 10"^^ 
9.0 X lO"^'* 


-1.5 
-3.0 


0.8 
1.5 
1.5 
1.5 


0.70 
0.30 
0.08 


0.306 

0.628 
0.275 
0.050 
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Fig. 1. — The dotted-line in this F\ vs Xrest 
plot represents the composite quasar spectrum of 
ZK97. The straight short-long dashed line rep- 
resents a power-law fit, in the longer wavelength 
region A > 1200 A {Xrest), of the continuum un- 
derlying the emission lines. The thick solid line 
overlaying the ZK97 data is our simulated com- 
posite Model A, which assumes a single power-law 
with oc to describe the intrinsic quasar 
SED and a tenuous H i absorption screen defined 
by eq. (1) with the parameters listed in Table 1. 
The short-dash line represents Model B, in which 
the Hi distribution corresponds to cq. (2) with 
zp = 0.7. The thin long-dashed line illustrates 
how deep the discontinuity at Xobs. = 1216 A is in 
models without z'q dependence (i.e. without cav- 
ity, the curve shown was obtained by setting /3 = 
in Model A). 




Fig. 2. — The dotted-line in this F\ vs Ar-est plot 
represents the composite spectrum of radio-quiet 
quasars by TZ02. The straight short-long dashed 
line represents a power-law fit (of index —0.72), in 
the longer wavelength region A > 1200 A {Xrest) of 
the continuum underlying the emission lines. The 
thick solid line overlaying the TZ02 data is our 
composite Model C, assuming a single power-law, 
Fi, oc for the intrinsic quasar SED, and a 

tenuous H i absorption screen described by eq. (3) 
with the parameters Hsted in Table 1 {zp = 0.3). 
The staircase thin continuous line represents the 
transmitted spectrum for an individual quasar of 
redshift zq = 1.0, assuming cq. (3) [Model C] 
and an ideal detector with sensitivity at all wave- 
lengths. The natural logarithm of the flux jump 
between the two filled squares is Tugo (sec defi- 
nition in §3.2.2). The thin long-dashed line illus- 
trates the discontinuity at Xrest — 1216 A if we set 
Zp to zero in Model C (i.e. without cavity). 
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Fig. 3. — The lower solid line is the calculated 
transmission curve for the quasar 0308—193 at 
zq = 2.756 as a function of the observer-frame 
wavelength for Model C [see eq. (3) and Table 
1]. The dotted lines shows the continuation of 
the transmission curve outside the HST-FOS spec- 
tral window used on that occasion (1315-2400 A). 
Starting from the right, the first, second and third 
ramps, seen in the transmission curve, are domi- 
nated by the opacity of Lya, Ly/3 and Ly7 lines, 
respectively. The line labeled "contribution..." il- 
lustrates the small contribution of photoionization 
to the total opacity. The upper solid line is the 
transmission curve for a lower redshift quasar at 
ZQ = 0.33. 
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Fig. 4. — Behavior of the depth of the flux dis- 
continuity, Tiieo, evaluated at 1160 A (Xobs.) (c-f- 
§ 3.2.2), as a function of quasar redshift zq for the 
Models A, B, C and D. 
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Fig. 5. — The dotted-line in this F\ vs Xrest plot 
represents the composite spectrum of radio-quiet 
quasars by TZ02. The straight short-long dashed 
line represents the intrinsic broken power-law SED 
assumed in Model D. It has a break at 1200 A 
at which point the power-law index changes from 
—0.72 to —1.57. The thick solid line overlaying the 
TZ02 data is Model D using the parameters listed 
in Table 1. As in Fig. 2, the thin continuous line 
represents the transmitted spectrum of an individ- 
ual quasar with zq = 1.0 (using Model D param- 
eters). The thin short-dashed line illustrates the 
discontinuity at Xrest = 1216 A if we set zp to 
zero in Model D (i.e. without cavity). 
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